Medium Modifications and Production of Charmonia at LHC 
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A previously constructed transport approach to calculate the evolution of quarkonium yields and 
spectra in heavy-ion collisions is applied to Pb-Pb(y / s=2.76 ATeV) collisions at the Large Hadron 
Collider (LHC). In this approach spectral properties of charmonia are constrained by euclidean 
correlators from thermal lattice QCD and subsequently implemented into a Boltzmann equation 
accounting for both suppression and regeneration reactions. Based on a fair description of SPS and 
RHIC data, we provide predictions for the centrality dependence of J/tj) yields at LHC. The main 
uncertainty is associated with the input charm cross section, in particular its hitherto unknown 
reduction due to shadowing in nuclear collisions. Incomplete charm-quark thermalization and non- 
equilibrium in charmonium chemistry entail a marked reduction of the regeneration yield compared 
to the statistical equilibrium limit. 

PACS numbers: 25.75.-q, 12.38.Mh, 14.40.Pq 



I. INTRODUCTION 

In recent years it has been realized that early ideas of 
associating charmonium suppression in ultrarclativistic 
heavy-ion collisions (URHICs) with the deconfinement 
transition ll I are less direct than originally hoped for, 
see Refs. 0tJ) for recent reviews. While there is little 
doubt that a sufficiently hot and dense medium created 
in URHICs induces a large suppression of primordially 
(or would-be) produced charmonia, coalescence processes 
of charm and anti-charm quarks in the later stages of the 
medium evolution are now believed to play an important 
role in a quantitative description of charmonium observ- 
ables. The regeneration reactions are, in fact, dictated by 
the principle of detailed balance and are operative once 
a charmonium bound state can be supported in the heat 
bath [5l4TTj]. The onset of regeneration for each char- 
momium state in the cooling process of the bulk matter 
is thus dictated by its "dissociation temperature", T^ lSB , 
which in general will be different for ground (^—i] c ,J/tp) 
and excited states (ty=ip' ,Xe, etc.). Furthermore, the in- 
elastic reaction rate, Ty, of each charmonium controls 
its approach to the equilibrium limit. Both reaction rate 
(inelastic width) and equilibrium limit depend on the 
temperature and thus change as the system evolves in 
a heavy-ion reaction. The balance between primordially 
produced and regenerated charmonia can therefore be ex- 
pected to be sensitive to the spectral properties of char- 
monia in the Quark-Gluon Plasma (QGP). 

In Ref. fl2| we have implemented basic properties of 
charmonium spectral functions (mass, width and binding 
energy) into a Boltzmann transport equation. The spec- 
tral properties have been estimated from an in-medium 
potential model [l3[ and checked against euclidean corre- 
lators from thermal lattice QCD. Since the definition of 
an in-medium two-body heavy-quark potential remains 
an open problem to date, we have investigated two limit- 



ing cases representing a "weak-binding scenario" (WBS) 
and a "strong-binding scenario" (SBS) which are char- 
acterized by a small (1.2 T c ) and large (2T C ) Jjvb dis- 
sociation temperature, respectively. The corresponding 
phenomenological analysis of J/ yj production in URHICs 
at SPS (Vs=17.3AGeV) and RHIC (yi=200 AGeV), 
solving a rate equation in a thermally expanding fire- 
ball background, revealed that both SBS and WBS can 
be compatible with currently available data, despite the 
factor of ~10 difference in collision energy. The key dif- 
ference is the composition of the final yield, especially in 
central A- A collisions, where the regeneration (primor- 
dial) component is much larger in the WBS (SBS) than 
in the SBS (WBS). However, since the individual exci- 
tation functions (-y/s dependence) of the primordial and 
regenerated components are quite opposite (decreasing 
vs. increasing, respectively) [14]], it is hoped that an- 
other factor of ~10 increase in collision energy is able to 
disentangle at least the two limiting scenarios. An initial 
investigation of this possibility is the main purpose of this 
paper, by providing predictions for inclusive J/vj produc- 
tion in Pb-Pb collisions at -^5=2760 AGeV as recently 
conducted at the Large Hadron Collider (LHC). Another 
important point is in how far chemical and thermal off- 
equilibrium effects cause a deviation from the equilibrium 
limit as predicted, e.g., in the statistical hadronization 
model 



,s predic 
15, 16]. 



The main uncertainty in our predictions will be asso- 
ciated with the input for the open-charm cross section 
(as well as its shadowing in Pb-Pb). Here, preliminary 
LHC data from p-p collisions will be utilized in connec- 
tion with theoretical estimates for shadowing effects. The 
charged particle multiplicity, dN c h/dy 1 has already been 
measured and will be used to constrain the total entropy 
in the system at each centrality. The total multiplicity 
is not only important for the initial temperature and its 
evolution but also determines the volume of the system 
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at given temperature which affects the charm-quark den- 
sity (for given total number of charm pairs, N c£ ) and thus 
the equilibrium limit of charmonia. 

In the following we hrst recapitulate our framework for 
calculating the evolution of charmonia in URHICs and 
the connection to their spectral properties (Sec. [IIJ) , dis- 
cuss the open- and hidden-charm input quantities for Pb- 
Pb at ^5=2760 AGeV (Sec. HI)), evaluate results for the 
centrality dependence of inclusive J/tp yields (Sec. HV|) 
and conclude (Sec. IV) . 

II. RATE EQUATION AND SPECTRAL 
PROPERTIES OF CHARMONIA 

In this section we briefly recall the main features 
of constructing the equilibrium properties of charmonia 
( i Sf=J/ijj, Xc an d ip') in the Quark-Gluon Plasma (QGP) 
and hadron gas (HG), and their im plem entation into a 
thermal rate equation following Ref. [12| (which contains 
a more detailed account of our approach). Upon inte- 
grating over the spatial and momentum dependence of 
the Boltzmann equation one arrives at a simplified (av- 
eraged) rate equation for the number, AT^(t), of state & 
in a heavy-ion collision, 

^ = -r*(T) [N 9 -N?(T)] . (1) 

The dissociation (and formation) rate, T^, is computed 
in the so-called quasifree approximation for the processes 
p + ty P+ c + c induced by light partons, p, of the heat 
bath (u, d, s anti-/quarks and gluons) pjj]. Besides the 
trivial T dependence induced by the number density of 
the light partons, the charmoniun binding energy, 

e B (T) = 2m* c (T)-m*(T) , (2) 

plays an important in that states with larger eb are char- 
acterized by a smaller r$(T); this ensures a realistic hi- 
erarchy between ground and excited states, and is re- 
sponsible for the main difference in T^,{T) between the 
SBS and WBS. The binding energies, £b(T), in these two 
scenarios are taken from T-matrix calculations 13] using 
1QCD results for either the internal ([/gg(r;T)) or free 
energy (Fqq(t;T)) as interaction potential, respectively. 
Within each scenario the in-medium charm-quark mass is 
fixed according to 2m* (T) = U q q(oo;T) or F q q(oo;T), 
which also determines the charmonium mass, m^(T), via 
Eq. The latter figures into the statistical equilibrium 
limit of each state, 

Nt\T) = j 2 c (N c - c ;T) V FB J -j0-f(p- T ) , (3) 

where d^ denotes the spin degeneracy, /* (p; T) the Bose 
distribution and Vfb the (time-dependent) fireball vol- 
ume. The charm-quark fugacity, 7 C (T") = lc{T), is de- 
termined by assuming a fixed number, N cS , of charm- 
quark pairs in the fireball for a given collision centrality, 



as following from primordial production in binary N-N 
collisions (the underlying cc cross section is one of the 
main input parameters discussed in Sec. IIIII below). For 
a given Vfb and temperature, j c is calculated from rel- 
ative chemical equilibrium among the open- and hidden- 
charm states in the system, where the former are either 
charm hadrons in the HG phase (using vacuum masses) 
or charm quarks with mass m*(T) in the QGP phase 
(charmonia have a negligible impact on j c at the tem- 
peratures of interest here). To account for incomplete 
thermalization of the charm-quark spectra in the QGP, 
we adopt a relaxation-time approximation to correct the 
statistical equilibrium limit of charmonia [l7| . The equi- 
librium limit figuring into the rate equation ([1} thus takes 
the final form 

N c ^ = K(t) < at , K(t) = 1 - exp(-r/T c e <i) . (4) 

The thermal relaxation time r^ q has been taken as a fit 
parameter which controls the regeneration contribution 
in the rate equation. The resulting values, r c cq ~ 2-4fm/c 
for the WBS and SBS, respectively, are on the same or- 
der as microscopic calculations of charm-quark transport 
which lead to a fair description of heavy-quark observ- 
ables at RHIC [l8|. The impact of the TZ- factor on the 
J/ip regeneration yield at SPS and RHIC is appreciable, 
and quite comparable to a schematic coalescence model 
where the limiting cases of thermalized and initial power- 
law spectra have been studied (l9| . Incomplete charm- 
quark thermalization will also play a significant role in 
our predictions for LHC as discussed below. 

To check the equilibrium properties of charmonia 
against 1QCD "data", we have "reconstructed" spectral 
functions by combining a relativistic Breit-Wigner ansatz 
for the bound-state part (using M^(T) and T^(T) as de- 
scribed above, as well as a polestrength factor Z-$(T) 
to mimic the wave-function overlap at the origin) with 
a nonperturbative continuum (with threshold 2m* (T)). 
Pertinent euclidean correlator ratios turn out to be rather 
stable with T, deviating from 1 by ca. ±10% for both SBS 
and WBS, roughly consistent with 1QCD. The vanishing 
of Z^(T) is used to quantify the dissociation tempera- 
ture, Tf iss , above which the gain term in the rate equa- 
tion is switched off. 

III. INPUT CROSS SECTIONS AT ^=2760 GEV 

In this section we collect the inputs required for the 
initial conditions of the rate equation in the charm and 
charmonium sector, as well as for the thermal fireball 
evolution. 

For the total open-charm cross section in p-p colli- 
sions at -y/s = 5.5 TeV recent calculations using next-to- 
lcading order (NLO) perturbative QCD (pQCD) obtain 
2.5-3.5mb [2(| using m c =1.5GeV, compared to 5-7. 5mb 
in earlier calculations with m c =1.2GeV [2l|. Prelimi- 
nary LHC data at 7 TeV [22| indicate a cross section 
close to the upper end of the pQCD predictions, which 
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is around 10 mb. Reducing this range by ca. 1/3-1/2 to 
extrapolate to y/s — 2.76 TeV, we estimate <Jp P {y/s = 
2.76TeV) ~5-8mb, in agreement with fixed-order-next- 
to-leading-log (FONLL) pQCD predictions [HHl]. Con- 
verting this to a midrapidity density 21 1 (which amounts 
to dividing by a factor of 7-8 [2j |2j, (24[), we arrive 

at ^f(\/s = 2.76 TeV) ~ 0.7-1 mb. Our calculations 
reported below are therefore conducted with a default 

value of ^(Vs = 2.76 TeV) = 0.75 mb (which leads to 
a reasonable (J/tp)/(cc) ratio when evaluating available 
J/tp cross sections, see below). The extrapolation to Pb- 
Pb collisions is done via standard binary-collision (N co \\) 
scaling using an inelastic p-p cross section of 65 mb ac- 
counting. At LHC energies the nuclear modification of 
parton distribution functions in the nucleon (shadowing) 
is predicted to suppress the charm cross section signifi- 
cantly. We will investigate this effect by performing cal- 
culations where, according to Ref. [2l|, the input charm 
cross section is reduced by up to 1/3 for central Pb-Pb, 
with a centrality dependence as estimated from Refs. [25j . 

For the charmonium cross section p-p measure- 
ments around midrapidity are available from CDF 
(Vs=1.96TeV) [il and ALICE (^5=7 TeV) Q with 

da J/i> 

= (3.4 ± 0.3)/x6 and (7.5 ± 2)fxb, respectively. In 
addition, charm measurements at lower energies suggest 
an approximately constant (energy-independent) frac- 

tion of J/tp to open-charm cross sections, J P / —f 2 - — 



dy 

dat'J 



(0.5-1) %, see, e.g., Ref. [15|. Thus we take — — 
2.76 TeV) = 4 /ib, which, with the above charm cross 
section, implies a J/ip-to-cc fraction of 0.53%. When in- 
cluding charmonium shadowing in Pb-Pb we assume a 
suppression by 1/3 in central collisions, consistent with 
recent estimates at y/s—5.5ATeV at midrapidity [27j us- 
ing the color-evaporation model; this is somewhat smaller 
than the factor of 0.5 suppression found in Ref. [28[ us- 
ing EKS98 29] or nDSg [301 shadowing, or a suppression 
by 2/3 or more predicted by the Color Glass Condensate 
approach HH, Cold- Nuclear-Matter (CNM) effects 
other than shadowing are neglected. Recent empirical 
evaluations of the nuclear absorption cross section in p- 
A collisions indicate a marked reduction with increasing 
\/s, suggesting it to be very small in the LHC energy 
regime. This is naively expected since the passage time 
of the nuclei at LHC has become very small. We also 
neglect the Cronin effect. An initial-state p t broaden- 
ing could be induced by shadowing effects which are ex- 
pected to be reduced with increasing p t , but we neglect 
this in the present work. It turns out that our results 
for the J/V> nuclear modification factor are not sensitive 
to CNM effects on the J/tp except for rather peripheral 
collision where the primordial component is sizable. 

For the thermal medium evolution, the main empiri- 
cal input is the total entropy, S, in the fireball which we 
assume to be conserved. We calculate the entropy den- 
sity using a hadron-resonance gas equation of state (EoS) 
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FIG. 1: (Color online) Time profile of temperature in central 
(6=0) Au-Au ( v / s=0.2ATeV, A^t ~ 380; dashed line) and 
Pb-Pb (Vs=2.76ATeV, N paTt ~ 400; solid line) collisions us- 
ing empirical charged-hadron multiplicities, dN^/drj ~ 750 
and 1750, respectively, as input. 



at chemical freezeout (s^ m ~ 6fm~ 3 with T c hem = 
180 MeV). The 3- volume is then adjusted to obtain a 
charged-hadron multiplicity of dN c h/drj ~ 1600 for 0- 
5% central Pb-Pb at ^=2.76 ATeV (N palt ^ 380) [33j], 
and the measured centrality dependence is employed for 
non-central collisions [34j |. Using a QGP quasiparticle 
EoS for T > T c = T chcm = 180 MeV, we connect the 
two phases via a first-order transition, which allows to 
construct the temperature evolution, T(t), based on a 
cylindrical isotropic fireball volume expansion, Vfb(t), 
using s(T) = S/Vfb(t). With a formation time of 
to = 0.2 fm/c the initial temperature amounts to To = 
610 MeV, compared to 330 MeV at RHIC. The lifetime 
of the QGP phase (where most of the J/tp chemistry is 
operative) is also significantly larger, tqgp — 6 fm/c vs. 
3 fm/c at RHIC, see Fig. HJ 



IV. CENTRALITY DEPENDENCE OF J/tp 
YIELDS IN PB-PB( v / 5=2760 GEV) 



With charmonium and open-charm input cross sec- 
tions, as well as the fireball evolution at each centrality 
fixed, we solve the rate equation (JTJ for the 3 charmonia 
VP = J/tp, ib' and \c and obtain the final inclusive J/tp 
yield using standard feeddown fractions, i.e., 32% and 
8% from ip' and Xc in pp collisions, respectively. In addi- 
tion we account for B-meson feeddown (Bfd) according 
to Tevatron data (26|, amounting to ca. 10% of the inclu- 
sive J / ip yield in pp (without any modifications of its pt 
dependence, i.e., 6-quark energy loss), see, e.g., Ref. [35j . 
We display the J / tp yields in terms of the nuclear modi- 
fication, either as a function of collision centrality (char- 
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acterized by the number of nucleon participants, -ZV par t), 

N Jfij {N palt ) 



RAA(N part ) 



JVcou(JVpart) N™^ ' 



(5) 



or at fixed centrality as a function of fireball evolution 
time, 



Raa(t) 



N con (N pa , Tt ) Nf fi> ' 
or transverse momentum, 



RAAipt) 



dNj/^/dpt 



N con (N paIt ) dNf^/dpt 



(6) 



(7) 



Let us first examine the time dependence of the ex- 
clusive J/ip yield as a function of fireball evolution time, 
displayed in Fig. [2] for central Pb-Pb in the weak- and 
strong-binding scenario (upper and lower panel, respec- 
tively), using = 0.75 mb without shadowing. In 
both scenarios the large initial temperature leads to a 
strong and rapid suppression of the primordially pro- 
duced J/ip's. The somewhat surprising feature is that 
the final regeneration yield is very similar in both sce- 
narios, despite the different binding energies and dissoci- 
ation rates at given temperature. The time dependence 
of the regenerated component furthermore shows that the 
cc coalescence occurs at rather different times or, equiva- 
lently, temperature, namely in the regime where the gain 
term in the rate equation first becomes operative, i.e., for 
Tdiss < 1-2(2)T C , corresponding to r > 4(l)fm/c in the 
WBS (SBS). It is precisely in this regime where the dis- 
sociation (and thus formation) rate is large, facilitating a 
rather rapid approach of the J/tp abundance toward its 
equilibrium limit. However, the "cooking" only lasts for 
a duration of ca. At ~ 2 — 3 fm/c in both scenarios, after 
which the reaction rate becomes small so that the yield 
stabilizes. 

The centrality dependence of the final J/ip yields (in- 
cluding feeddown) in Pb-Pb collisions is displayed in 
Fig. [3l again for both WBS (upper panel) and SBS (lower 
panel) and for the same conditions as in the time depen- 
dence shown in Fig. [2] (i.e, iV C oii-scaled charm cross sec- 

tion of -jP' = 0.75 mb without shadowing), except that 
bottom feeddown (Bfd) and formation-time effects (fte) 
are now included as detailed in Ref. [35[. One finds that 
the similarity of the two scenarios not only holds for cen- 
tral but also for semicentral and even peripheral colli- 
sions. For iVp ar t ~ 200 the composition of the total J/ip 
yield in terms of primordial and regenerated components 
is quite comparable to central Au-Au collisions at RHIC, 
as is the charged-particle multiplicity. However, due to 
a smaller initial overlap volume at LHC (at identical 
dN c h/dy) the initial temperature is significantly larger 
than at RHIC. In addition, charm production is larger 
and thus regeneration effects are more pronounced. 
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FIG. 2: (Color online) Time dependence of the J/ip nu- 
clear modification factor in central (b=0) Pb-Pb colli- 
sions at v / s=2.76 ATeV (participant number TVpart — 400) at 
1/5=2.76 ATeV within the WBS (top) and SBS (bottom). No 
shadowing effects are included; only the exclusive J/ip contri- 
bution is included in the numerator of Raa (t) , Eq. ([6} , while 
its denominator includes the 40% feeddown fraction, to fa- 
cilitate the comparsion with subsequent plots (consequently, 
Raa(t=0)=0.6). Bottom feeddown and formation-time ef- 
fects are not included. 



A striking feature of our predictions is a rather large 
deviation of the final J /ip yield from the values pre- 
dicted by the statistical model of hadron production (up- 
per dash-dotted lines in Fig. [3]). By definition, these 
abundances refer to an equilibrated hadron gas at hadro- 
chemical freezeout (which is naturally identified with the 
hadronization transition) , with a fixed charm-quark num- 
ber enforced by the fugacity, 7 C , in Eq. ([3]). The equilib- 
rium limit, A^, q , in our rate equation ([T]) basically coin- 
cides with the statistical equilibrium limit, iV|, tat , since 
the thermal off-equilibrium 7\L-factor has become close to 
one at the end of the mixed phase (recall Fig. [2]). While 
the observed light- and strange-hadron abundances in 
central AA collisions are in excellent agreement with the 
statistical model [36J, our calculations suggest that the 
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FIG. 3: (Color online) Centrality dependence of the nuclear 
modification factor, Eq.(JSJ), for the inclusive J/ip yield (in- 
cluding feeddown) in Pb-Pb(y / s=2.76 ATeV) collisions within 
weak- and strong-binding scenarios (upper and lower panel, 
respectively). 



inelastic charmonium reaction rates are too small in the 
vicinity of T c , especially on the hadronic side, to establish 
relative chemical equilibrium in the charm/onium sector. 
Consequently, the J/ip yields following from the kinetic 
approach are well below the statistical-model limit (by 
about a factor of ~3), resulting in a suppression below 
one in the nuclear modification factor. 

As is well known, in the grand canonical limit the equi- 
librium ^ number depends quadratically on the number 
of open-charm pairs in the system, oc 7^ oc AT?-. 
Therefore, regeneration (or statistical production) are 
quite sensitive to the open-charm cross section. In Fig. 0] 
we display our results for R^(N paT t) in the SBS in- 
cluding a (pt -independent) shadowing correction to both 
N cS (b) = N con (b) N*£ and Ny(b), as described in Sec. [Oil 
As expected, the reduction of the N cS by 1/3 in central 
Pb-Pb suppresses the regeneration J/ip yield by a fac- 
tor of ca. (2/3) 2 ~ 1/2, so that the total Raa decreases 
from ^0.8 to 0.45 (no shadowing is applied to the bot- 
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FIG. 4: (Color online) Nuclear modification factor for 
the inclusive J/ip yield (including feeddown) in Pb- 
Pb(y / s=2.76 ATeV) collisions within the strong-binding sce- 
nario including a shadowing correction of the iV co ii-scaled 
open-charm and J/ip cross sections by a factor of 2/3 for 
iV part ~400. 



torn feeddown contribution). This underlines again the 
importance of an accurate open-charm input cross sec- 
tion for in- medium charmonium physics at LHC. 

Next, we investigate the transverse- momentum spec- 
tra of inclusive J/ip production, displayed in Fig. [5] 
for peripheral and central Pb-Pb collisions (upper and 
lower panel, respectively). For simplicity we restrict 
ourselves to predictions without shadowing corrections. 
At both centralities the regeneration yield generates a 
marked enhancement of the spectra at low p t . The 
maximum structure is more pronounced in central col- 
lisions due to larger regeneration and larger suppres- 
sion of the primordial production; the crossing between 
the two components occurs at pt — 5-6 GeV, compared 
to pt ~ GeV for peripheral collisions. The latter is 
reminiscent of central Au-Au collisions in the SBS at 
RHIC 12], although there no significant maximum de- 
velops at low p t (the charm cross section is too small and 
the suppression of the primordial component is stronger 
in central Au-Au(ys=0.2 ATeV) than in peripheral Pb- 
Pb (y/s—2. 76 ATeV). We note, however, that shadowing 
corrections are expected to be most pronounced at low 
Pt (i.e., small values of x — 2p t /y/s at midrapidity) and 
fade away at high pt- Since shadowing suppresses both 
regeneration and primordial production, our calculations 
in Fig. [5] may be viewed as an upper limit for the maxi- 
mum at low p t . 

Finally, we perform a calculation for the so-called 
central-to-peripheral ratio of J/ip mesons, which is simi- 
lar to Raa but normalized to a peripheral centrality bin, 



R C P (A^part) 



A^(AWt)/AUi(AWt) 

/yperi /Ty-peri 
ly ih / Jv coll 



(8) 



A first measurement of this quantity in Pb-Pb at LHC 
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FIG. 5: (Color online) Transverse-momentum dependence of 
the nuclear modification factor, Eq.lJT]), for inclusive J ftp's 
(including feeddown) in peripheral (upper panel) and cen- 
tral (lower panel) Pb-Pb( ^=2.76 ATeV) collisions within the 
strong-binding scenario (no shadowing corrections accounted 
for). 



has recently been reported by the ATLAS collabora- 
tion [37[ where the peripheral bin corresponds to a 40- 
80% centrality selection. In addition, only muons with 
Pt > 3 GeV were accepted, implying that ca. 80% of 
the reconstructed J/ip carry a transverse momentum of 
Pt > 6.5 GeV. We mimic these conditions in our cal- 
culation by representing the 40-80% centrality bin with 
an average iV par t — 40 and using a momentum cut of 
p t > 6.5 GeV for all J/^'s (neglecting shadowing). The 
latter implies that the regeneration contribution is es- 
sentially gone, recall Fig. [3] so that our calculated Rcp 
only involves primordially (suppressed) J/^'s, including 
bottom feeddown and formation-time effects in the QGP 
suppression reactions (35| . The trend of our postdiction 
(without any adjustment of parameters) is fairly com- 
patible with the ATLAS measurement, cf. Fig. [6] We 
emphasize again that the high-pt cut evades conclusions 
about regeneration mechanism other than that they are 
not operative in this regime. 
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FIG. 6: (Color online) Central-to-peripheral ratio, Rcp, 
for inclusive J/ip yields (including feeddown) in Pb- 
Pb(y / s=2.76 ATeV) collisions within the strong- and weak- 
binding scenario (solid and dashed line, respectively; no shad- 
owing included). A p t > 6.5 GeV cut has been applied to all 
J/ip's in the calculations to compare to ATLAS data [3?| . 



V. CONCLUSION 

Employing a kinetic rate-equation approach which in- 
corporates in-medium effects extracted from charmonium 
spectral functions we have computed inclusive J/tp pro- 
duction in Pb-Pb collisions as recently conducted at 
the LHC. Our approach reproduces existing charmo- 
nium data at SPS and RHIC energies where the two 
main model parameters had been adjusted (strong cou- 
pling constant in the dissociation rate and thermal off- 
equilibrium correction in the equilibrium limit). Our per- 
tinent predictions for LHC indicate a large suppression of 
primordial charmonia (by a factor of ~10 in central Pb- 
Pb) and thus a predominance of the regeneration yield for 
participant numbers V par t > 100 — 150. However, the fi- 
nal yield stays well below the predictions of the statistical 
model, by about a factor of ^3, implying an R J Jj{ below 
one even for an input charm cross section of 0.75 mb in 
pp collisions; shadowing corrections on this quantity fur- 
ther suppress the final J/ip yield. Somewhat surprisingly, 
both strong- and weak-binding scenarios give very similar 
results, which we attribute to the rather limited time win- 
dow (At ~ 2 — 3 fm/c) during which the gain term in the 
rate equation is active (below Tdi SS ) and the inelastic re- 
action rate is large. In the strong-binding scenario J/^'s 
form at significantly higher temperatures and thus at ear- 
lier times in the fireball evolution than for weak binding. 
This may lead to discernible differences in the finally ob- 
served elliptic flow, which, however, will have to be mea- 
sured with rather good precision. We also performed cal- 
culations of R J Jp in comparison to very recent ATLAS 
data, finding a suppression roughly compatible with ex- 
periment. However, due to the large-pt cut in these data 
(pt > 6.5 GeV), our result is not sensitive to regeneration 
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contributions which only figure at lower p t < 6(3) GeV 
for central (peripheral) collisions. Upcoming LHC data 
are eagerly awaited to test our predictions and give fur- 
ther insights into charmonium production mechanisms in 
hot/dense QCD matter. 
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